Spin selective transport through hehcal molecular systems 
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Highly spin selective transport of electrons through a helically shaped electrostatic potential is 
demonstrated in the frame of a minimal model approach. The effect is significant even for weak 
spin-orbit coupling. Two main factors determine the selectivity, an unconventional Rashba-like spin- 
orbit interaction, reflecting the helical symmetry of the system, and a weakly dispersive electronic 
band of the helical system. The weak electronic coupling, associated with the small dispersion, 
leads to a low mobility of the charges in the system and allows even weak spin-orbit interactions 
to be effective. The results are expected to be generic for chiral molecular systems displaying low 
spin-orbit coupling and low conductivity. 
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Introduction— The concept of spintronic devices oper- 
ating without a magnetic field has been proposed some 
time ago for solid state devices in which the spin-orbit 
couphng (SOC) is large [J, 0. Recently, a new type of 
magnet-less spin selective transmission effect has been 
reported [^-Q. It was found that electron transmission 
through chiral molecules is highly spin selective at room 
temperature. These findings are surprising since carbon- 
based molecules have typically a small SOC that can- 
not support significant splitting between the spin states, 
splitting which is thought to be essential for any spin de- 
pendent property. Although it has been found both in 
theory 0-01 and experiments that there is a cooper- 
ative contribution to the value of the SOC, so that this 
quantity may be larger in molecules or nanotubes than 
in a single carbon atom, the values calculated or exper- 
imentally found are still relatively small [2"11|, e.g. few 
meV for carbon nanotubes [lo| . Hence, even including 
this cooperative contribution, the values obtained for the 
spin polarization (SP) in electron transmission through 
chiral molecules [6| seem to be too high and cannot be 
rationalized by such SOC values. 

Recently, a theoretical model based on the first Born 
approximation in scattering theory has been proposed for 
explaining the spin selectivity of chiral molecules 13| . Al- 
though the results are in qualitative agreement with the 
experimental observations, they could not explain them 
using reasonable SOC values. 

In what follows, a model is presented to describe elec- 
tron transmission through a helical electrostatic poten- 
tial (see Fig. 1). Although the model does not claim to 
fully catch the complexity of the experimentally studied 
DNA-based systems @, [ij] , it highlights the role of some 
crucial parameters, which can determine the experimen- 
tally observed high SP. The key factors in the model that 
allow for the high spin selectivity are: i) Lack of parity 
symmetry due to the chiral symmetry of the scattering 



potential; ii) Narrow electronic band widths in the helical 
system, i.e. the interaction between the units composing 
the helical structure through which the electron propa- 
gates is relatively weak. Moreover, a physically meaning- 
ful estimation of the SOC can be obtained by taking into 
account that first, like in the solid state, in the present 
study the electric field acting on the electron needs to in- 
clude the effective infiuence of all the electrons belonging 
to a molecular unit 13, 15 [ , and second, due to proximity 
effects, the Coulomb interaction between the transmitted 
electron and the atoms in the molecular unit can scale as 
1/R for short distances R. 

Model and Methodology— We consider the Schrodinger 
equation for a particle moving in a helical electro- 
static field. Analytical results for such fields have 
been derived in Ref. [l^- For the sake of simplicity, 
approximate expressions valid near the z-axis will be 
used (only x and y components will be considered. 
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FIG. 1. A charge q in spin state a is moving along through he- 
lical electric field. The parameters a, b and Az are the radius 
and the pitch of the helix and the spacing of the z-component 
of the position vector of the charges distributed along it, re- 
spectively. The helical field Ehdix induces a magnetic field B 
in the rest frame of the charge and hence influences its spin 
state. 
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the z component only contributes when considering the 
fuh three-dimensional problem, see below ): Ehciix — 
-EQj2ijgi,j{z){cos{QjAz),sm{QjAz)). Reve, gij{z) = 
(1 + jAz)/a]2)-3/2 and Q = 27r/6 with b being 

the helix pitch and a the helix radius, see Fig. [T] The 
index m — 0, - ■ ■ , Mq — 1 runs along one helical turn and 
labels the z-coordinate of the Mq molecular units placed 
along the helix. The index n = — Lo/2, ••• ,Lo/2 {Lq 
being the number of helical turns) connects sites which 
differ in their z-coordinate by b [16|. We note that the 
considered helical potential is assumed to be related to 
the charge distribution along the stack of molecular units 
building the helical structure; hence the factor Eq is pro- 
portional to the local charge density. 

For a charge moving with momentum p through the 
helix, the field EhcUx induces a magnetic field in the 
charge's rest frame, from which a SOC arises: Hso = 
Aa(p X Ehciix)- The SOC strength is A = efi/Arri^c^ and 
<T is a vector whose components are the Pauli matrices 
CTj;, cTj^, fjz- The general problem is three-dimensional; 
however, in order to get first insights into the behavior 
of the SP, we will assume Px — Py — 0,Pz 7^ 0, so that 
the Schrodinger equation takes the form [17]: 



f{z) 

r{z) 



C 

X{z) = Ex{z) 
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Here, x{z) — {x' [z) , i^))'^ is a spinor, ^(z) = — 



AQjAz 



g^,Jiz), f{z) = dz'i'iz), and U{z) the 
helical electrostatic potential. The terms ~ f{z),f*{z) 
are introduced to make the Hamiltonian hermitian in 
the continuum representation. The SOC parameter a = 
HXEq (with dimensions of energy x length) depends on 
the effective charge density through Eq. The problem 
posed by Eq.[T]can be written as an effective two-channel 
nearest- neighbor tight-binding model ^17i] : 

N N-1 
^ = ^ Uncl^aCn,a + ^ ^ i^n,aC,i+l,a + h.C.) 

N 

+ J2 {cl^Wn,,nC,n,i+C+^^W^,„Cn,t)+Hleads- (2) 



The operators {Cn,a, ct^a}n=i,. 



,N,, 



create or destroy. 



respectively, an excitation at the tight-binding site n with 
spin index a. The only non-zero elements of the inter- 
channel coupling matrix W are given by [l7| : Wn,n = 
-a/(nAz), Wn,n+i = Q*(nAz)/2Az, and W„+i,„ = 
— + l)Az)/2Az. Further, the matrix satis- 
fies W^,^ = -(Wn,m)* for n ^ m, and W^„^„ = (W„,„)*. 
The hopping V should in general be estimated on the 
basis of a first-principle calculation of the electronic cou- 
pling for a given system. However, we will consider it 
as a free parameter, whose order of magnitude for helical 



organic systems is expected to lie in the range of few tens 
of meV (e.g. for DNA, electronic structure calculations 
yield values of the order of 20 — 40 meV [l^)- Finally, 
the operator Hieads includes the semi-infinite chains to 



the left (L) and right (R) of the SO active region [17 [. 
A schematic representation of this two-channel model is 
shown on the top panel of Fig. [2l 

Transport properties— We focus on the spin-dependent 
transmission probability, T[E), of the model Hamilto- 
nian given by Eq. [21 as a function of the electron's in- 
jection energy E. The problem can be considered as a 
scattering problem where a finite-size region (with non- 
vanishing SOC) is coupled to two independent L (left)- 
and two independent i?(right)-electrodes, each electrode 
standing for a spin channel and being represented by a 
semi-infinite chain, see Fig. [5J T{E) encodes the influ- 
ence of multiple scatte ring events in the SOC region; us- 
ing Landauer's theory [20| we obtain 17i]: 



r(£;) = r^(r^|Git,^tl 



+ ^i\GlUN-\\'^) 



+ rf|G 



(3) 



^) spin- flip scattering, see 
the trans- 



In Eq. [31 G na,mv{E) with (T, v — 'tii- ^rc matrix elements 
of the retarded Green's function of the SOC region in- 
cluding the influence of the L- and i?-electrodes. The 
individual contributions in Eq. [3] can be related to dif- 
ferent transport processes without (e.g. r:^r:|^|Gi-|-^Ar-|-p) 
or with (e.g. r:^rf |Gi^ 
Fig. [H Notice that tup{E) and tdown{E) 
missions for the up and down channels respectively, as 
defined by Eq. [3]-, contain contributions arising both 
from direct transmission without spin-flip as well as 
spin-flip. An energy-resolved SP for different initial 
spinor states can be defined as: P{E) = {tup{E) — 
tdown{E))/T{E). The energy-average SP {P{E))e = 
P{{tup{E)),{tdown{E)),{T{E))) will also be used. We 
focus only on electron-like contributions {E < 0) and 
on energies \E\ > ksT « 23 meV, so that {■■■)£ = 

j:^fdE{...). 

Results— A crucial parameter in the model is the SOC 
coupling a. Realistic values are obviously very diffi- 
cult to obtain 2l|, [22] , since a is not simply the atomic 
SOC, but contains the influence of the charge distri- 
bution in the system via the field factor Eq. For the 
sake of reference, a rough value of Eq for DNA may 
be estimated along the following lines. A single DNA 
base is considered as composed of discrete point-like 
charge centers A, representing the atoms. We associate 
with each center A at position Ra a Gaussian-shaped 
charge distribution of width w ^ 0.3 — 0.4 nm and with 
strength given by an estimated atomic charge density 
po for C, N, and O atoms (considered as spheres with 
a radius of the order of the corresponding covalent ra- 
dius). The local field of this charge distribution, Eq = 
-(l/47reo)(9/9r) / d^r' p{r' - RA)\r - r'\~^ , canhe com- 
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FIG. 2. Top panel: Schematic representation of the tight- 
binding model, see Eq. [2] The two channels interact via the 
SOC (framed region). To the left and right of the spin scatter- 
ing region, both channels are independent and are modeled 
by semi-infinite chains. Bottom panel: Energy dependence 
of the SP P{E) for Lq—S helical turns, and for injected elec- 
trons polarized with their spin pointing up (Pio), down (Poi), 
or unpolarized (-Pii). A spin-filter effect takes place only for 
energies near the band edges, where all SPs have the same 
sign. Notice also that near the band edges the SP has oppo- 
site signs for electrons {E < 0) and holes {E > 0), though 
P{E) is not exactly antisymmetric. Parameters: a — 5 meV 
nm, V = 30 meV, f/o = 3 meV. 



puted analytically [l^ and it scales for R = j?- — <C w 
like Eo ~ {NopQ/ATTto){w/2^)'^R~'^ (Eq has been mul- 
tiplied by a factor A^o 10, the number of atoms in a 
base, to approximately account for other charge centers. 
For R/w ~ 0.3 - 0.4, values of a = hXEa ~ 1.87 - 2.35 
meV nm can be obtained. In the calculations, a ^ 2 — 6 
meV nm have been used. Though the previous discussion 
provides a very rough estimate, it highlights the need of 
considering the influence of many charges through pq and 
Nq as well as proximity effects (short-distance scaling of 
Eq) in the estimation of a. 

Fig- HI presents the energy dependent SP for different 
incoming spin states when the spin is pointing up (10), 
down (01) or the electrons are unpolarized (11). The 
coupling a was assumed to be 5 meV nm. Although 
this value is larger than the previously estimated one, it 
serves to illustrate the behavior of the model in a clear 
way. In the case of (10) and (01) states, the interesting 
energy windows are those where both SPs have the same 
sign, which indicates that the outgoing state will always 
have the same SP independently of the initial condition. 
This behavior occurs mainly for energies near the band 
edges. A similar situation is found for the (11) state, see 
Fig. Near the band center, Pio{E) and Poi{E) have 
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FIG. 3. Different components of the transmission 
tupiE),tdown{E) and T{E) as defined in Eq. [3] and for the 
same parameters of Fig. [2] Focusing on electron-like contri- 
butions, it is only near the lower band edge {E < —22 meV) 
where a positive SP for all incoming states (10), (01), and 
(11) is obtained, see also Fig. [21 



opposite signs and hence the SP depends on the incoming 
spin state. The average SPs, as defined above, amount 
to approximately (Pio)£; = (^oi)e = {Pii)e w62%. 

Figure |3] shows the corresponding spin- resolved trans- 
missions. Notice first, that the states (10) and (01) corre- 
spond to cases where one of the incoming spin channels is 
decoupled from the system by setting the corresponding 
= 17[, and hence, the total transmission cannot 
exceed one. For (11) both channels are open and the 



maximum transmission is 2. 

In the top panel of Fig. |31 for (10) and (01), we find 
some degree of spin-dependent back-scattering, which 
is reflected in the different total transmissions T{E) 
for each polarization. In what follows, for the sake 
of reference, only the behavior in the energy window 
[-2V,-kBT],kBT - 23 meV is discussed. For the (10) 
state, transmission without spin flip is dominant in this 
energy region, and this leads to the positive SP. How- 
ever, for (01), spin-flip processes become dominant in the 
same energy region, and hence the outgoing up-channel 
acquires a larger weight. As a result the SP for (01) 
is also positive. This behavior is closely related to the 
chiral symmetry, which basically manifests in the spe- 
cial structure of the W,W^ matrices. For the (11) state, 
bottom panel of Fig. |3l the outgoing up-channel clearly 
dominates the transmission in the considered energy win- 
dow , thus indicating that for unpolarized electrons back- 
scattering and spin-flip of the down-component will ul- 
timately lead to a positive SP. A similar analysis can 
be performed for the hole- like energy region E > 0. In 
general terms, SP may occur either by spin-flip (with no 
net change of the total transmission) or by spin selec- 
tive back-scattering. The results of Fig. [3] suggest that 
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FIG. 4. 2D plot of the energy average SP {P{E))e as a func- 
tion of both the hopping parameter V and the SOC a. Only 
for small V a relative large SP is found. With increasing elec- 
tronic coupling, larger SOC strengths are required to get a 
sizeable SP. 



influence the SP. The present study indicates that CISS 
may be a very general phenomenon, existing in chiral sys- 
tems having low SOC and low conductivity, and hence 
may play a role in charge transport through biosystems. 
The effect could also be of great interest to control the 
spin injection efficiency in the context of semiconductor- 
based spintronics by interfacing chiral molecules with 
semiconductor materials. 
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both processes are playing a role; their relative contribu- 
tion to the SP turns out however to sensitively depend 
on the specific energy window considered. The selectiv- 
ity found in this model relates to two special features of 
the chiral system: (i) the symmetry of the field which 
translates into an unconventional SOC, and (ii) the nar- 
row electronic band width in chiral organic systems. The 
term band width serves only as a keyword for the aver- 
aged value of the coupling matrix elements, V, between 
neighboring molecular states mediating charge motion. 
As shown in Fig. |4l the size of the hopping parameter 
strongly affects the energy average SP, ultimately lead- 
ing to {P{E))e for large V. For small hopping, 
however, the SP can achieve very large values by only 
a moderate increase of the SOC a. The interplay be- 
tween a and V seems related to the relatively long time 
(roughly proportional to hV^^) the electron will spend 
in the conducting channel in a real system, allowing for 
the SOC to become more effective. 

Conclusions— The present study based on a generic 
model sheds new light on a chirality-induced spin selec- 
tivity (CISS) effect. It suggests that beyond the symme- 
try itself, CISS depends on the organic molecules being 
poor conductors. Weak electronic coupling along the he- 
lical structure is expected to lead to low mobility of the 
electrons through the system and allows enough time for 
the SOC, although being weak, to influence spin trans- 
port. The effect depends on the electron momentum 
and once the electrons have kinetic energy above kgT, 
the SP increases and becomes weakly energy dependent. 
One open issue for further inquiry is the influence of the 
electrode-molecule interface. If the electrodes are mag- 
netic, spin-dependent tunnel barriers emerge, which may 
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